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Abstract 

The difference between induced by box diagram quantity e and experimentally 
measured value of e is determined and used to obtain the value of e with high 
precision. Present day knowledge of CKM matrix elements (including B-factory 
data), allows us to obtain from the Standard Model expression for e the value of 
parameter Bk'- Bk = 0.89 ± 0.16. It turns out to be very close to the result of 
vacuum insertion, Bk = 1. 
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1 Introduction. 



It is well known that CP - violation in K° — K° mixing is described by the parameter e. 
Within the SM, this parameter is given by box diagrams. It depends in particular on the 
CKM matrix elements, to which vertices of box diagrams are proportional. On the other 
hand, the experimentally measured parameters are e and e'. e and e' enter the measured 
ratios of decay amplitudes of kaons into mi states. These amplitudes are superpositions 
of amplitudes A(K° — > (mv)j) = A I e tSl of kaon decays into states with definite isospin 
7 = 0,2, Aj are weak amplitudes, 5i are strong rescattering phases of 7r-mesons. The 
parameter e can be expressed as [T]: 



e = e + % - - . (1) 

Within the SM and in the standard parametrization of CKM matrix, ImAo originates 
from the so-called strong penguin diagrams. Amplitude A 2 also has an imaginary part 
which originates from electro- weak penguin diagrams. That is why ImA » ImA 2 . 

Taking into account that the phases of e and e are approximately | pQ, from Eq.Q 
we obtain: 

, , 1 ImA n (2) 



y/2 ReA c 



ImA 



The estimation of — — — was done in [2]. This term appears to be a 5 — 9% correction 
ReAo 

to the value of e in Eq.(J2J). Provided that we have estimated the right-hand side of Eq.(j2*|) 
with the help of Eq. (jl3|) we can determine parameter Bk, which parameterizes hadronic 
matrix element. Of course, for this purpose we need to know the values of CKM matrix 
elements that enter Eq. ffTS"]) . 

The parameters p and fj of CKM matrix appear to be constrained without using 
the value of e in the fit. Thus we perform the fit of CKM matrix parameters without 
using constraint from e in it. Then we determine from Eqs.(j2*j). (J13|) . Our result is 
B K = 0.89 ±0.16. 

This result is close to the result of vacuum insertion: Bk = 1. As discussed in 
the insertions of 7r-mesons states should be taken into account. These insertions form a 
sign-alternating series, who's terms depend on the cutoff momentum of 7r-mesons. This 
cutoff can be reasonably chosen to be 200 — 500 MeV (at larger virtualities 7r-mesons do 
not exist). Then the sign-alternating series converges quickly, and one can take only first 
two terms. Thus taking into account the insertions of 7r-mesons states lowers B^, and 
the agreement with our result improves further. 

The lattice result of Bk calculation is Bk = 0.87 ± 0.06 ± 0.14 quenc h jl]. We see that 
our result is very close to it. 

The paper is organized as follows: in Section |21 we discuss various estimations of the 

value of -— — ■ . In Section El we perform the fit of CKM matrix parameters without using 

ReAo 

constraint from e in it. In Section 0] we determine Bk and compare it with other results 
of calculation of Bk- Finally, we make our conclusion in Section |SJ 
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2 Estimation of the numerical value of e — e. 

In this section we review the estimation of e — e |2]. We discuss the following three 

methods. First, one can use the experimental data on CP-violation in semileptonic Kl- 

decays, namely parameter 5 L . This method possesses large uncertainty and also at the 

e' 

level of two sigmas contradicts the experimental value of — . Second, one can obtain the 

e 

e' 

lower bound on e — e from the experimental value of — . This lower bound is important 

e 

in understanding the relative magnitude of the second term in Eq.(J2J), it turns out to 

be > 5%. Third, we use the results of direct computation of ImA in the ratio — — — , 

ReA 

substituting the experimental value of ReA . This gives us a reliable estimate of — — 

ReA 

with moderate error,which we use in the bulk of the paper. 

First, we estimate the value of e from the experimental results on CP- violation in 
semileptonic K L decays: 

f T(K L -> 1+V7T-) - T(K L -> 1-97T+) ^ 2Rei 



T(K L -> l+v-K-) + T(K L -> Z-P7T+) 

Sl 
2 cos 6 ' 



(3) 



where <fi = arg(e). 

Now let us substitute the experimental data. For we use (ft = (43.50 ± 0.05)° jSJ. 
World average value of Sl, published in jHj, contains new KTev result: Sl = (3.307 ± 
0.063) x 10~ 3 . Na48 collaboration recently obtained: S L = (3.317 ± 0.100) x 10" 3 [Zj. 
Averaging these two numbers we get: Sl = (3.310 ± 0.053) x 10 -3 . This leads to the 
following value of e: 

|e| = (2.282 ± 0.037) x 10" 3 . (4) 

From Eq.© with the help of Eq.© we can find the corresponding value of -^: 

ReA 

— — - = 0.03 ± 0.56 x 10" 4 . 5 
ReA 

We will show below, that this number almost contradicts the present experimental 
value of - = (1.67 ± 0.26) x 10" 3 0. 

Second method of estimation of -, which gives the lower bound on it, uses the 

ReA 

e' e 1 
experimental value of — . The expression for — is usually presented as follows ^: 



1 = _ c i(&2-5 ) 1 

e y/2 e 



ImA 2 ImA 

w 

ReA ReA 



(6) 



Let us neglect the term proportional to I171A2 in Eq. fBj) . which comes from the EW 
penguins. Taking into account that (So — S<i) exv = 42 ± 4° [8 a , we obtain the following 

expression for from Eq.©: 

ReAn 
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e' 



ImA _ _V2\e[e^_ 
ReA w e 

Substituting experimental values from 0, we get: 

1.67 ±0.26) x KT 3 , w = 0.045, |e| = 2.284(14) x 10~ 3 =>• 

1 ^ = -(1.2±0.2)x 10 - (8) 

In this way we get the following value of |e|: 

|e| = 2.37(2) x 10" 3 . (9) 

Since, according to Eq.flnj), the contribution of EW penguins partially cancels that of 

QCD penguin, the value I— — -^1 = (1.2 ± 0.2) x 10~ 4 should be considered as a lower 

ReA Q 

bound on I— — -^1 and Eq.Q is a lower bound on lei. Thus the central value of — — 
ReA ReA 

obtained from semileptonic -K^-decays, Eq.fl^J), almost contradicts the experimental value 

off: 

e J A 

Finally, the reliable way to estimate — , result of which we will use in the next 

ReA 

sections is to use the experimental value of ReA and theoretical value for ImA . 

Calculation of ReAo and IitlAq, as well as ReA2 and ImA2, has a long history. The 
calculation of ReAo and ReA2 was performed in order to explain the AI = ~ rule in kaon 
decays and the calculation of ImAo and ImA?, - in order to explain the observed value of 

e 

In this paper we perform the calculation of ImA to the following accuracy: the 
Wilson coefficient is calculated to LO and hadronic matrix element is calculated in naive 
factorization approximation. The details are presented in Appendix, and here we only 
quote the result: 



We note that the results of computation of ImA (see jHj - [Hj and refs. therein) 
performed by a large number of people lie in the same ballpark. 
Finally, from Eq. fjlO)) we can determine the value of e: 

|e| = (2.51 ±0.07) x 10~ 3 . (11) 
This number is our final result, and we will use it in Section 



3 Fit of the parameters of CKM matrix 

We use in our fit of the CKM matrix experimentally measured values of modulus of matrix 
elements V^V^V^V^V^s, and also sin2a, sin2/3, sin2"f and Ams d . Note that we 
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do not use e in fit, since we plan to determine the value of Bk with the help of the fit 
results. 

We assume these experimentally measured data to be normally distributed. Also the 
theoretical uncertainties are treated as normally distributed. Let us note that other people 
treat theoretical uncertainties in other way ^B] , [TTj . 

The table of input parameters looks like: 



Parameter 


Value 


Standard Deviation 




v ud 


u 


0.9738 


0.0005 




V us 





0.2200 


0.0026 




v ub 


U 


0.00367 


0.00047 




v cd 




0.224 


0.012 




v cs 




0.996 


0.013 




v cb 




0.0413 


0.0015 


sin2a 


m 


-0.21 


0.46 


sin2j3 |S] 


0.736 


0.049 


sm27 [11] 


0.69 


0.58 



The x expression which we minimize looks 



X 2 (A,X,p,r]) 



ike: 



+ 



+ 



■\rtheo \r 

v cs v _t 

a V cs 

' sin2(3 theo 



exp 



( sin2a 




theo 



a v cd 
sin2a> exp \ 



" sinloL 



sin2f3 



\ ' sinlfi J \ ~ siril'y ) 

where theoretical expressions depend on four Wolfenstein parameters: A, A, p and f\. 
Expression (|T2*|) was minimized varying them. 
Here are our results: 



+ 



(12) 



a 

fj[ de 9\ 
[deg] 
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= 100 ±5 
23 ±2 
57 ±5 



A = 0.224 ± 0.002 ™ [deg] 

A = 0.82 ±0.03 
p = 0.22 ±0.04 
fj = 0.34 ±0.02 

X 2 /n.d.o.f. = 8.1/5 . 

For comparison, we present the results of the fit, made by CKMfitter Group ^H] and 
UTfit Collaboration \TT\: 





CKMfitter 


UTfit 


A 
A 

P 
fj 


0.226 ±0.002 
0.80l°;°l 

0.19185? 

0.36l°;°I 


0.226 ±0.002 

0.17 ±0.05 
0.35 ±0.03 





CKMfitter 


UTfit 


a [deg\ 
R\deg\ 


wis 

Z6.6_2 .0 

6211° 


94 ±8 
23.2 ± 1.4 
61.6 ±7 
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4 The value of B K 



From the results of the fit, presented above, we can extract the value of Bk- For this 
purpose we use the theoretical expression for |e|, first obtained in [20] • It has the following 
form: 



,^ eo| = B K ( Vcc miIm[(V cs V^] + Vu miI(£)Im[(V ts V td 



G 2 F m K f K 



l2V2ir 2 Am K 



,m 2 w . 



+2r) ct m 2 \n(^)Im[V cs V; d V ts V; d }). (13) 

Here /({) = { ^±± _ Jf±U { = M . Quark mas8e8 are m „ = L2 ± 

0.2 GeV [5], m t = 178.0±4.3 GeV [21], m w = 80.42±0.04 GeV 0. The QCD corrections 
were calculated to leading order in (201 : Vcc = 0.6, r\ tt = 0.6, = 0.4. The next- 
to-leading order calculation changes slightly r\ tt and r\ ct and changes considerably r\ cc : 
r] cc = 1.32 ± 0.32 [22], Vtt = 0.574 ± 0.01 [2H], r/ c t = 0.47 ± 0.04 |2U - The kaon decay 
constant extracted from the K + — > ji + u decay width equals: fx = 160.4 ±1.9 MeV jS]. 
The K L - K s mass difference is Am K = (3.483 ± 0.006) x 10~ 15 GeV 0. Fermi constant 
Gp = 1.16639(1) x lO- 5 GeV- 2 0. 

Now we equate this expression to the value of |e| = (2.51 ± 0.07) x 10~ 3 from Eq.(|llj). 
substituting all experimental numbers and the results of the fit. This leads to the following 
value of Bk '- 

B K = 0.89 ±0.16 (14) 
Note that it is close to the result of vacuum insertion: Bk = 1- 



5 Conclusions 

We have extracted the value of Bk using the fitted values of CKM matrix elements and 
the estimated difference between e and e. Our result is Bk = 0.89 ± 0.16. It appears to 
be close to the result of vacuum insertion, Bk = 1, while lattice result is simply the same: 
B K = 0.87 ± 0.06 ±V.U quench H]. 
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A Estimation of the value of IitiAq from QCD pen- 
guin diagram. 

Let's estimate using experimental value of ReA and evaluating the value of ImA . 

The latter will be evaluated to the following accuracy: the LO Wilson coefficients will be 
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used, and hadronic matrix element will be calculated in naive factorization approximation. 

As it is well known transitions with AS = 1 are due to the 4-quark effective Hamilto- 
nian, for the first time derived in [25j: 

6 

H AS =x = V2G F sin 9 C cos 9 C £ cA (15) 

i=l 

The so-called penguin operator O5 dominates in amplitudes K° — > (7T7r)/ = o 

5 = s L 7 At A a rf L (M JJ 7 / ,A a M H + d R ^X a d R ) (16) 

Below we present the detailed derivation of the coefficient function C5 in one loop 
approximation. 




W W 



Figure 1: Diagrams, which contribute to the penguin operator. 

It is convenient to perform calculation in the unitary gauge. Each of the three diagrams 
(see Fig. ^) is infinite, however the sum appears to be finite. We will use the dimensional 
regularization (d = 4 — 2e), in order to regularize divergent integrals: 



9 2 9l 


1 


J=U,C,t 


f d d k 9fj.u 


kfj,ki/ 


2 


q 2 


(27l) d k 2 




9 2 9 2 s 


1 


J=U,C,t 


f d d k 9^u 


kfikv 
M 'w 


2 




(2vr) d k 2 




9 2 9l 


1 


J=U,C,t 


f d d k 9/j.v 


M 'w 


2 


q 2 


(2vr) d k 2 





J*' 

""J \ / 

' " ' ' ,_ " '■ 1 qi — rrij p2 — rrij 2 \ 2 j 

M<2 = U / 7^—2 — ^^Tr^w^WLd (vvt^Vv! V jd v* s , 

M = M + M 1 + M 2 . (17) 
The effective Hamiltonian is equal to: 

Hpenguin iM. (18) 

The result of the calculation can be presented in the following form: 
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H penguin = - £ ^S L G p {q) yrf L (VV7 P yVvj V^*, 

G p = Gi7 P + -Jy- (G 2 p2l P Pi + Gzpipip + G A p 2 p 2p + G b pxp2 P + G 6 p 2 p lp ) . (19) 



Dimensionless functions (i = 1, 2, 3, 4, 5, 6) depend on the values Xj = jp-, q , m 2 . 

m 2 ,, where rrij is the up-quark mass in the loop. We suppose the external s and d quarks 
to be on mass shell. 

Let us neglect d-quark mass. In this approximation the non-zero contribution into 
operator H penguin is given by the terms with formfactors G\, G^ and Gq. It is convenient 
to introduce new variables P = p\ + p 2 and q = p\ — p 2 . 

s L G p —d L = G x s L -i p —d L + ^ -s R —d L . (20) 

Here the term proportional to q p will not contribute to operator H penguin , since q p x 

(Vv7p4fVv) = ( m r — m r) (^p2^" l Pi) — 0. The quantity P p s R ^d L should be expressed 
through the magnetic formfactor with the help of the following equation: 

sa^qJ^-dL = ^qusilvlv - lulp) y^L = ^(^s^yc^ ~ P^R^-di). (21) 
With the help of Eq. (j2Tfl from Eq.® we obtain: 

tj sr^ 9 2 9s 1 / t - , . . h - , \ 7 A a 

H pe nguin = - 2^ J^2~ nSLlp—d L + ijr^rn s q v s R a ^—dL ip r ^ p —ip r V jd V js) 



(22) 



where fi and f 2 are equal to: 



2 

mi 



fl — G\ + - - (Gq + G4J, 



2M W 



h = (23) 

It is sufficient to calculate the formfactors G4 and Gq in the zero order in q 2 , m 2 s . How- 
ever, the formfactor G\, as it follows from last equations, should be calculated, including 
terms proportional to m 2 , q 2 . From equations (JUj), calculating appropriate integrals, we 
get: 

G\ = -RlTTo 1" -^2 



M 2 W M^' 
_ 7x 4 + 14x 3 - 63x 2 + 38x + 4 + 6(16x - 9x 2 - 4) \nx 
1 ~ 24(1 - xY : 

-5x 4 + 14x 3 - 39x 2 + 38x - 8 + 18x 2 In x 
2 ~ 8(1 -xY ' 



s 



2x 4 - 14x 3 + 45x 2 - 38x + 5 + 6(1 - Ax) \nx 
G4 = 



6(1 -x) 4 

llx 4 - 14x 3 + 27x 2 - 38x + 14 + 6(8a; - 9x 2 - 2) \nx 

12(1 -x) 



G 6 = TT^ 3i " " • ( 24 ) 



Substituting these formulas into equations (|23|) . we obtain: 

_ 7x 4 + Ux 3 - 63x 2 + 38x + 4 + 6(16x - 9x 2 - 4) lux q 2 
h ~ 24(1 -x) 4 M^' 

5x 4 - 14x 3 + 39x 2 - 38a; + 8 - 18a; 2 lux 



Finally, let us rewrite equation (j22|) in the following way: 

Hpenguin = ~ a o 2*1 2 f F ^Llu^-d L + lF 2 m s ^SrCT^ ^d L J ^rT/^^VVj'd^* 

i=«,c,t 487r iW W V ^ 1 Z / * 

7x 4 + Ux 3 - 63s 2 + 38x + 4 + Q(lQx - 9x 2 - 4) Ins 



24(1 -x) 4 



5a; 4 - 14a; 3 + 39a; 2 - 38x + 8 - 18a; 2 Ins 
F2 = h = 8(1^? • (26) 

As the admixture of gluons in K and 7r mesons is small, the contribution of magnetic 
moment operator in (j2Bj) is negligible [2*5] . 

Substituting m c = 1.2 GeV, m t = 178.0 GeV, M w = 80.42 GeV for the formfactor F 1 
we obtain: 

Fi(x « 1) « — lnx + -, 

6 

Fx (x c ) ^8.58, 
F x (x t ) « 0.550, 

*i (oo) = ^ « 0.292. (27) 
Formula ()26j) can be rewritten with good accuracy as: 

U „ fan ■ a a ( Qs i m c . JmV cd V* s a s 

Hpenguin ~ V 2Gjr sin C cos C - — In — + z— — — — In — - 5 , (28) 

\ 12n fi 2 ReV c dV c * 12tt m z c J 

where instead of m u the characteristic hadronic scale fi (this time "low" normalization 
point) is substituted. 

Thus the real and imaginary parts of C5 are equal to: 

In order to understand at which virtuality a s should be taken in these expressions 
leading logarithms should be summed up. This was done for the real part of coefficient 
function in the paper 
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Rec 5 = (x?- 48 (-0-039 X °2 S + 0.033 X ° 2 A2 + 0.003 X J - 12 + 0.003 x 2 - a3 ) + 
+xr°' 24 (-0.014 X f - 0.001 x^ 42 - 0.014 X2 ai2 + 0.029 x 2 ~ a3 )) , (30) 



while for imaginary part in paper the following result was obtained: 

j mC5 = /m ^f (0.0494 x?- 85 " 0-0280 x? 42 + 0.0116 xf 013 - 0.0330 x^ 35 ) * 
x (0.8509 x°- 8 + 0.0091 x^ 42 + 0.1222 xj " 12 + 0.0178 x 2 °' 3 ) , (31) 

where X i = i! r ! H, X2 = 

A - 1 a s (m w )' AZ a s (m c ) 

Numerical analysis shows that with a good accuracy expression for Rec$ can be written 

as: 



127T 8 >2 



ite* = -^log(^). (32) 



On the other hand Imc$ at the scale // at which a s (n) = 1 has the value 

/mC5= ife^ x013 - (33) 

The expression for ImAo can be written as: 



ImAo = V2G F sm8 c cos6 c Im(c 5 ) < (7rn) I=0 \O 5 \K > (34) 

In order to get the value of ImA we must calculate hadronic matrix element of 
penguin operator. It was evaluated in the framework of naive quark model in see 
also (23: 



4:V6 m 2 K mlU ( Ik 



m 



< (7rn) I=0 \O 5 \K° >= ^ TTJ n ^L (1 + ^) _ ! (35) 



9 m s (m u + m d ) \ f n m\ , 

Substituting experimental numbers Gf = 1.166 x 10~ 5 GeV, sin 9c = 0.22, cos 9c = 
0.95, m n = !35MeV,m K = A97MeVJ w = 130MeVJ K = 160MeV,m a = 700MeV and 
quark masses m s = 130MeV,m u = 3MeV,m d = 7MeV, we get: 

ImA = -11 x KT^GeV. (36) 
Finally, dividing it by experimentally measured ReA = 3.33 x 10~ 7 GeV, we obtain: 

ImAn , 

— -£ = -3.2 x 10- 4 . 37 

ReAn 
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In order to estimate the theoretical error for ImA we propose the following method: 
to take two values of /i, corresponding to a s (/i) = | and a s (jj) — |, to calculate at 
each /z and from these boundary values get ± error for . 

Via the proposed method we get our final result: 

^ = -< 3 - 2 -> x 10 " 4 < 38 > 

This number is rather stable with respect to the variation of /i. Our result confirms 

statement maid in |2Zj: QCD penguin results in the value of — in ballpark of the experi- 
mental data. 

On the other hand the real part is very sensitive to /i, as can be seen from Eq.()32|). 
The expression for ReA has the following form: 



ReA = V2G F sin 6 C cos 6 c Re(c 5 ) < (nn) I=0 \O 5 \K° > . (39) 

Substituting numbers and again taking fi at which = | and a s (/x) = | we get: 

ReA = (1.2^Qg) x 10~ 7 GeV. The central number is approximately 3 times smaller than 
the experimental , but theoretical uncertainty is large. 
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